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Abstract 

The North American Linear Colhder Detector group has developed a simulation program 
package based on the ROOT system. The package consists of Fast simulation, the recon- 
struction of the Full simulated data, and physics analysis utilities. 
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1 Introduction 



For the various studies of the future Linear Colhder experiments, the North American Linear 
Colhder Detector group (LCD) has constructed a simulation facihty. It consists of two kinds 
of detector simulators, Full and Fast simulators, designed for detailed detector studies and 
physics analyses, respectively. We use the GISMO package [1] for the Full simulation. Using 
the digitized outputs from the Full simulation, we reconstruct charged tracks and clusters. 
The Fast simulation is based on parametrized position and energy smearing, and it makes 
tracks and clusters directly from the generated particle information. 

In the LCD group, we have developed two object-oriented program packages based on 
JAS/Java[2] and Root/C++[3], for the reconstruction of the Full simulated events, the Fast 
simulation and the event analysis. In this report, we introduce the simulation and analysis 
package based on Root/C++. 

Fig.l shows the LCD simulation and analysis flow using ROOT. We use a standard 
HEPEVT format using the FNAL StdHep 4.06 I/O package [4], for the generator output and 
the detector simulator input. The output format of the Fast simulator is the Root binary, it 
is also possible to be linked with the event analysis directly. The output format of the Full 
simulator is the SIO binary and is converted to the Root binary. Then its reconstruction 
and the event analysis are done based on the LCD Root system. 
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Figure 1: The LCD Simulation/Event analysis flow. Boxed processes are based on the ROOT 
system. Several event generators are also ROOT-based. 
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Table 1: Detector parameters for LCD Small and Large detectors. 





Small 


Large 


Vertex Detector 


CCD 


CCD 


Impact parameter resolution 


3.O//08.7///psin^/^ 


2.5/1® 8.7 iJ,/p sin^/^ 


Central Tracking 


Si micro strips 


TPC 


Momentum resolution (High) 
(Low) 


d/P'^ ~ 6 X 10-5 
5P/P - 0.4% 


S/P^ ~ 5 X 10-5 
SP/P ~ 0.1% 


Electromagnetic Calorimeter 
Barrel Inner Radius 


W/Si pads 
75 cm 


Pb / scintillator 
200 cm 


Endcap Inner Z 


150 cm 


300 cm 


Energy resolution 


5EjE ~ 12%/y/E + 1% 


SE/E ~ 15%/v^ + 1% 


Granularity 


20 mrad 


20 mrad 


Hadron Calorimeter 
Barrel Inner Radius 


Cu / scintillator 
140 cm 


Pb / scintillator 
250 cm 


Endcap Inner Z 


186 cm 


350 cm 


Energy resolution 
Granularity 


8E/E ~ 50%//E + 2% 
60 mrad 


SE/E ~ 40%/^/E + 2% 
60 mrad 


Coil Magnet 






Magnetic field 


6 Tesla 


3 Tesla 


Inner Radius 


100 cm (outside EM Cal) 


376 cm (outside HAD Cal) 



2 Simulation and Event Reconstruction 

There are two detector models presently under study by the LCD group, a Small and a 
Large detector. The detector parameters are hsted in Table 1. The Small detector is a 
compact detector featuring a very strong magnetic field and sohd-state tracking device. The 
calorimeter is also compact and made with Si-W. The Large detector has a much larger 
tracker (TPC) with lower magnetic field and Pb-scintillator calorimeter. Both detectors are 
designed to have a CCD Vertex Detector. These detector parameters are used for both Full 
and Fast simulations. 

In the Fast simulation, charged particles within the magnet field follow helical trajectories, 
and their momenta and positions arc smeared. Here the charged tracks are expressed by 5 
parameters, and are smeared with a 5x5 covariant error matrix. 

Electrons, photons, and hadrons produce clusters in the electromagnetic (EM) and 
hadronic (HAD) calorimeters. Here, one cluster is made from one particle. Energies and 
positions of clusters are smeared. We assume transverse position resolutions of 1cm/ y/E 
(electrons and photons) or 5cm/ \/E (hadrons). To consider the detector granularity and 
cluster width, which is typically a few units of Moliere radius, we merge the clusters when 
the angular separation between clusters is less than O^ax, where 9max is a size of the detector 
granularity. 

The position of the interaction position is also smeared. We assumed ax—(Ty—2 /im and 
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(7^=6 /xm. The determination of these values are described in Ref.[6]. 

In the reconstruction of the Full simulated data, we postpone the track reconstruction. 
Instead, we make charged tracks by smearing, using exactly the same procedure as in the 
Fast simulation, but we also apply a minimum tracker- hit cut. Calorimeter clusters are made 
by gathering the hits which are from the same particle. Energy and position of the cluster 
is obtained from the energy sum, and the energy- weighted average of associated Calorimeter 
hits, respectively. 

Comparing the Full and Fast simulations, the most significant difference is in the Calorime- 
try. It is important to improve the parameterization of the Fast simulator Calorimetry to 
have a more realistic detector response, which we hope to implement in the near future. 

3 Event Analysis Tools 

For the physics analysis, there are several useful tools in the LCD Root program. We 
provide Thrust finding and 3 kinds of Jet finding (based on JADE, JADE-E and DURHAM 
algorithms) programs. 

A topological vertex finding algorithm is developed in the SLD experiment [5]. Here the 
secondary vertices are reconstructed with charged tracks by searching the space points where 
track density functions overlap in the 3D space. The original SLD prepmort program, called 
ZVTOP, is translated into C++ and several parameters are set for the LCD studies. Based 
on the topological vertex finding, we get excellent performances on heavy-quark (both b and 
c) tagging and charge identification of the quark. The details are described in Ref.[6]. There 
are physics studies which use the topological vertex finding in these proceedings [7]. 

In the energy fiow analysis, neutral clusters are selected by the absence of a track and 
cluster association. For this analysis, we provide methods which extrapolate the particles to 
the cluster cylindrical radius. 

As a graphical tool, there is an event display for the LCD Root program. An example of 
the event view for a e'^e~ — > tt event is shown in Fig. 2 

4 Summary 

In this report, we briefly introduce the simulation and analysis tools based on ROOT for the 
LCD group. Although the tools are constantly being improved, they can be obtained via 
the URL, 

http : / /www-sldnt . slac . Stanford . edu/nld/New/Docs/LCD_Root/root . htm . 

Feedback from the users is highly welcomed. 
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Figure 2: An event display of a simulated 500 GeV e'^e — > event. Here the Fast simulator 
assuming the Large detector parameters is used. 
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